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Abstract 

Constraint of prompt photon data on the polarized gluon distribution is discussed in 
terms of uncertainty estimation for polarized parton distribution functions (PDFs). 
By comparing uncertainty of the double spin asymmetry Al r with expected statis- 
tical errors at RHIC, we found that the gluon distribution is effectively constrained 
j — , in the region 0.04 < x T < 0.2 with the data at transverse momentum p T = 10 — 20 

CN | GeV for center-of-mass energies ^fs = 200 and 500 GeV. 
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1 Introduction 



By the recent global analyses with the experimental data for polarized deep 
inelastic scattering (DIS), the polarized quark and antiquark distributions are 
determined well [1,2,3,4]. These distributions are obtained with enough ac- 
curacy to indicate that the quark spin content is smaller than prediction of 
naive quark model; AS = 0.1 ~ 0.3 1). These polarized parton distri- 
bution functions (PDFs) reproduce well the experimental data; however, the 
polarized gluon distribution Ag(x) cannot be constrained because of indirect 
and small contribution through Q 2 evolution and higher order correction at 
next-to-leading order (NLO). Furthermore, PDF uncertainty estimation indi- 
cated large uncertainty of the gluon distribution. It implies difficulty of the 
Ag(x) determination with only the polarized DIS data. 
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As a probe for the polarized gluon distribution, prompt photons will be mea- 
sured by the longitudinally polarized proton-proton collider at RHIC [5] . The 
gluon distribution contributes directly in the quark-gluon compton process 
(qg jq) at leading order (LO), and the process dominates in the whole p T 
region. The future asymmetry data contain useful information for clarifying 
the gluon contribution to the nucleon spin. Therefore, we are interested in 
the influence of the prompt photon data on the Ag(x) determination by the 
polarized PDF analysis. 

In this paper, we consider constraint of prompt photon data on the polar- 
ized gluon distribution. For evaluating the data constraint, we will compare 
the uncertainty of the spin asymmetry with the expected statistical error by 
the RHIC experiments. The asymmetry uncertainty coming from the polar- 
ized PDFs is estimated by the Hessian method, and it is comparable with the 
measurement error. In this comparison, the statistical error of the spin asym- 
metry plays a role of constriction for the Ag uncertainty via the asymmetry. 
In practice, it is indicated that uncertainties of the polarized PDFs can be re- 
duced by including new precise data for the polarized DIS in the Asymmetry 
Analysis Collaboration (AAC) [1]. Therefore, the same thing is expected by 
including the future data for prompt photon production. 



2 Uncertainty of the spin asymmetry 



The spin asymmetry A\ L is defined as a ratio of the polarized and unpolarized 
cross sections: A[ L = Aa 1 jo 1 . By the factorization theorem, the polarized 
cross section Act^^paPb — > lX) as a function of the transverse momentum p T 
is expressed by 

= J2 I dr] [ dx a [ dx b (1) 

dp T ~^ Jrj-bin J J 

x Af^(x a ,n F ) Af b B {x b ,fi F ) 
dAaJ b r 

X dp d ^ {x a ,X b , y/S,p T ,Tj, fl R , fl F ) , 

where Af a (x) is the polarized PDF of the parton a. We choose the AAC03 
PDF set [1]. Aalb is the partonic cross section (a + b — > 7 + X). In order to 
reduce theoretical uncertainty from the scale dependence of the cross section, 
the NLO corrections are taken into account. The NLO partonic cross sections 
for prompt photon production are completely known [6]. The renormalization 
and factorization scales are chosen the scale fj, F = fi R — p T . In addition, the 
cross section is integrated over the rapidity bin \q\ < 0.35, which corresponds 
to the acceptance of the PHENIX detector. The unpolarized cross section a 7 
is similarly computed with unpolarized PDFs and partonic cross sections [6]. 
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We choose the GRV98 unpolarized PDF set [7], which is also used in the 
AAC03 analysis. These cross sections are numerically calculated at center-of- 
mass (cm.) energy y/s = 200 and 500 GeV, respectively. 

In this study, the contribution from fragmentation is neglected. The contri- 
bution is associated with the collinear process for a scattered parton into a 
photon, and it can be diminished by using an isolation cut on the measured 
photon [8]. And so we should consider the isolation cut in this analysis [9]. 
For examining an effect on the polarized PDF uncertainties, we calculate the 
cross section for inclusive direct photon production process. 

The asymmetry uncertainty is obtained from uncertainty of the polarized cross 
section: 5A\ L = 5Aa J /a"'. The uncertainty SAa 1 comes from the polarized 
PDF uncertainties, and can be estimated by the Hessian method. The method 
based on eigenvectors of the diagonalized Hessian matrix is developed by 
CTEQ collaboration [10], and it is applied to estimate uncertainties of un- 
polarized PDFs [11]. We used the basic method, which is generally used by 
the polarized PDF analyses [1,2,3]. The uncertainty is given by 

where are optimized parameters in the polarized PDFs. Hij is the Hessian 
matrix which has information on the parameter errors and correlation between 
each parameter. The gradient terms of the cross section dAa 1 {p T )/dai are 
obtained as follows: 



dAa> 



ddi a b Jbin 



drj dx a dx a (3) 



x 



dAaJ b r 

X dp {x a ,X b ,y/s,p T ,Tj,ll R ,fl F ) . 

The gradient terms of the polarized PDFs can be derived analytically at initial 
scale, and these terms are numerically evolved to arbitrary scale Q 2 (= p 2 T ) by 
the DGLAP equation. Furthermore, the value of A% 2 determines a confidence 
level of the uncertainty. It is obtained so that the level corresponds to ler 
error of the normal distribution [1]. The uncertainty therefore can be directly 
compared with the statistical error of experimental data. 

For comparison with the asymmetry uncertainty, the statistical error of the 
spin asymmetry is estimated by 

P*\JC int O 
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where P is the beam polarization, Ci nt is the integrated luminosity, and a is 
the unpolarized cross section integrated over the p T bin. In this study, the a 
is computed by the bin size of 5 GeV interval, and it is taken the same bin 
size for both cm. energies y/s = 200 and 500 GeV. Furthermore, other values 
are chosed design values at RHIC [5]: P = 0.7 and Ci nt = 320 (800) pb -1 for 
y/s = 200 (500) GeV. 



3 Constraint of prompt photon data on Ag(x) 



First, we discuss predicted spin asymmetry and its uncertainty at x/s = 200 
GeV. In Fig. 1, the spin asymmetry by the AAC03 PDF set is compared 
to those by polarized PDF sets of BB (ISET=3) [2], GRSV01 (standard 
scenario) [4] , and LSS (MS scheme) [3]. These analyses used almost the same 
experimental data sets for the polarized DIS, and they obtained good agree- 
ments with the data. However, there are significant differences of the gluon 
distributions among them. 1 These differences are obviously reflected in varia- 
tions of the predicted asymmetries. The prompt photon process is sensitive to 
the behavior of the gluon distribution. Moreover, the asymmetry uncertainty 
is indicated in the same figure. Dotted curves show the uncertainty which 
comes from the polarized PDF uncertainties obtained by the AAC analysis 
with the polarized DIS data. We find that these predicted asymmetries are 
within the large uncertainty. These variations are caused by weak constraint 
of the polarized DIS data on the gluon distribution. The prompt photon data 
therefore are required for reducing this asymmetry uncertainty. 

In order to evaluate the gluon contribution to the asymmetry uncertainty, 



1 See, for example, Ref. [1,2] 
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Fig. 1. Comparison of the predicted spin asymmetries by different polarized PDF 
sets; AAC03(NLO), BB (ISET=3), GRSV01 (standard scenario), and LSS (MS 
scheme). The dotted curves show the asymmetry uncertainty from the PDF uncer- 
tainties of the AAC03 set. 
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we compute the asymmetry uncertainty excluding the Ag(x) uncertainty by 
assuming dAg(x)/da g = in eq. (3). Figure 2 shows the asymmetry uncer- 
tainties for \/s = 200 GeV. The solid curves show the current uncertainty 
by the AAC analysis, and the dotted curves do the asymmetry uncertainty 
except the Ag(x) uncertainty The significant reduction of the uncertainty in- 
dicates that the Ag(x) uncertainty is the dominant contribution to the current 
uncertainty. 

In addition, the asymmetry uncertainty is compared to the expected statistical 
errors at RHIC. The current uncertainty is much larger than these statistical 
errors. If these data are included in the global analysis, the uncertainty is 
roughly reduced to these errors. This suggests that the Ag(x) uncertainty is 
mainly improved. The prompt photon data therefore have strong constraint 
on the gluon distribution. 

On the other hand, the data constraint on the quark and antiquark distribu- 
tions is very weak. This is because that the asymmetry uncertainty without 
the Ag(x) uncertainty, which is composed of the Aq(x) and Aq(x) uncertain- 
ties, is significantly less than the statistical errors. Therefore, the data with 
such errors do not directly affect improvements of these uncertainties. 

However, as far as the antiquark is concerned, the uncertainty can be indirectly 
reduced because the antiquark distribution is strongly correlated with the 
gluon distribution in the global analysis. In practice, reduction of the antiquark 
uncertainty via the error correlation is indicated by the analysis with the fixed 
Ag{x) = at initial scale [1]. This fact suggests that the constraint on the 
gluon distribution indirectly affects the Aq(x) determination. In particular, 
it is not neglected when we perform flavor decomposition of the antiquark 
distributions. 



Next, we estimate a constraint factor for the Ag(x) uncertainty. By multiply- 
ing the gradient terms for the gluon dAg(x)/da g by the factor in eq. (3), the 
constricted uncertainty of the asymmetry is defined. From comparison of the 
asymmetry uncertainty with the expected statistical errors for \/s = 200 in 
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Fig. 2. Comparison of the asymmetry uncertainty SAJ^ L with the expected statistical 
errors for y/s = 200 GeV. The dashed curves show the asymmetry uncertainty except 
the gluon uncertainty: 5Ag(x) = 0. 
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the region 10 < p T < 20 GeV, the obtained factor is 1/18. 



In this study, the factors for the Aq(x) and Aq(x) uncertainties are neglected. 
This is simply because that these uncertainties are not directly constricted 
by these data. Moreover, the correlation effect on the Aq(x) uncertainty is 
not taken into account. The effect cannot be evaluated without including ex- 
perimental data in the global analysis. Since the Aq(x) contribution to the 
asymmetry uncertainty is already small, the uncertainty will be slightly mod- 
ified in this process. 

In Fig. 3, the constricted asymmetry uncertainties are compared with the 
expected statistical errors for yfs = 200 and 500 GeV, respectively. The solid 
curves show the asymmetry uncertainties which are obtained from the Ag(x) 
uncertainty multiplied by the constraint factor, and involve the Aq(x) and 
Aq(x) uncertainties. If x T {= 2p T /y/s) can be approximated by the Bjorken 
x(— x a = Xb) around central rapidity region, the data for yfs = 200 GeV in 
the region 10 < p T < 20 GeV constrain the gluon distribution in the range 
0.1 < x < 0.2. Although this comparison is in ideal condition that these data 
are put on the predicted asymmetry, this fact agrees with the results of the 
trial analysis including pseudo-data for Al L in Ref. [12]. These data are useful 
in determining the gluon distribution. 

In the region p T > 20 GeV, these data have rather weak constraint since 
the unpolarized cross section rapidly decreases as p T increases. The statistical 
errors depend on the p T bin size computing a in eq. (4). We should be careful 
about taking the bin size for the high p T data in order to constrain equally 
the gluon distribution over a wide x region. 

For the comparison at \fs = 500 GeV, we find similar behavior. In the region 
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Fig. 3. Comparison of the constricted asymmetry uncertainties with the expected 
statistical errors for y/s = 200 and 500 GeV, respectively. The solid curves show the 
constricted asymmetry uncertainty with the factor 1/18 for the Ag(x) uncertainty. 
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10 < p T < 20 GeV, the asymmetry uncertainty roughly corresponds to the 
statistical errors. This indicates that these data have the same constraint as 
those for y/s = 200 GeV, and constrain the gluon distribution in the range 
0.04 < x < 0.08. Above the region, the statistical errors are larger than 
the asymmetry uncertainty. It is noteworthy to mention here that the data 
constraint for yfs = 500 GeV is weaker than that for y/s = 200 GeV in the 
region 10 < p T < 20 GeV in spite of covering the same x T region. 

The reason for the weak constraint is that the unpolarized cross section for 
y/s = 500 GeV is less than that for y/s = 200 GeV in the same x T region, and 
the integrated luminosity is still insufficient to provide the enough constraint 
at high p T . Figure 4 shows the comparison of unpolarized cross sections for 
both cm. energies. These cross sections are indicated as a function of x T . In 
the region x T > 0.1, the cross section for y/s = 500 GeV is below one for 
200 GeV, and indicates similar behavior. In order to obtain equal constraint 
at the same x T , we need more luminosity than the design value at y/s = 500 
GeV. Therefore, the experimental data for y/s = 500 GeV are required as 
constraint on the gluon distribution at low x. The medium- x behavior should 
be determined by using the data for 200 GeV. 

Finally, let us turn to the Ag(x) uncertainty from the prompt photon data at 
RHIC. Figure 5 shows the polarized gluon distribution and its uncertainties 
at p T = 10 GeV. The solid curves show the Ag(x) uncertainty from the po- 
larized DIS data, and the shaded area shows the constricted uncertainty by 
comparison with the expected statistical errors. The uncertainty estimated by 
the constraint factor is reliable in the range 0.04 < x < 0.2. As a practical esti- 
mation, the uncertainty broadens gradually in the low- and medium- ;c regions 
where data do not exist. 2 Moreover, we note that increasing the asymmetry 
uncertainty with p T is due to uncertainty of the ratio of the polarized and 
unpolarized gluon distributions: 8Ag(x)/g(x). The uncertainty significantly 
increases with x due to lack of precise data for polarized DIS at large x. 

2 See, for example, Fig. 5b in Ref. [12]. 
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Fig. 4. The unpolarized cross sections for y/s = 200 and 500 GeV. 
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Fig. 5. The polarized gluon distribution and its uncertainty at p T = 10 GeV. Shaded 
area is the constricted gluon uncertainty. 

The Ag(x) uncertainty becomes the same order of magnitude as the quark 
uncertainty from the DIS data. By including future asymmetry data for the 
prompt photon process, the gluon distribution can be obtained with sufficient 
accuracy. 

In this study, polarization of the gluon distribution is not discussed. The cur- 
rent uncertainty in Fig. 5 indicates that we cannot rule out the possibility 
of zero or negative polarization. As another probe for the gluon distribution, 
the double spin asymmetry for ir° production has recently been reported by 
the PHENIX collaboration [13] . Since gg — > tt°X subprocess dominates at low 
p T , the cross section depends on (Ag) 2 ; therefore, the process is not sensitive 
to the sign of the gluon polarization. On the other hand, the prompt photon 
production which the qg compton process dominates is sensitive to the sign 
in the whole p T region. The gluon polarization is obviously reflected in the 
spin asymmetry. In this sense, the role of prompt photon data is of prime 
importance for determination of the gluon polarization. 



4 Summary 

In this paper, we have considered the uncertainty of the polarized gluon distri- 
bution for prompt photon production at RHIC. The uncertainty of the double 
spin asymmetry is estimated by the Hessian method. The asymmetry uncer- 
tainty mostly comes from the Ag(x) uncertainty. The large uncertainty implies 
the weak constraint of the polarized DIS data on the gluon distribution. By 
comparison with the expected statistical errors at RHIC, we indicate that the 
prompt photon data have the strong constraint on it. Furthermore, we suggest 
that the prompt photon data in the region 10 < p T < 20 GeV effectively con- 
strain the gluon distribution. The data of both cm. energies constrain it in 
the following regions: 0.04 < x < 0.08 at y/s = 500 GeV, and 0.1 < x < 0.2 at 
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200 GeV. For clarifying the gluon contribution Ag(= Jq dxAg(x)), the data 
covering a wide range of x are required. These experiments therefore play an 
important role in the Ag(x) determination. 
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